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Abstract. We study the properties of ultracompact jets 
in several prominent radio sources (Cygnus A, 3C 309.1, 
3C345, 3C395, 4C 39.25, and 1038+528 A), using the fre- 
quency dependence of observed position of the optically 
thick jet core. Frequency dependent offsets of the core po- 
sitions are used for calculating the luminosities, magnetic 
fields, and geometrical properties of the jets. Opacity effects 
are studied, for both synchrotron self-absorption and ex- 
ternal free-free absorption. Pressure and density gradients 
in the jets and in surrounding ambient medium are shown 
to be a primary factor determining the observed properties 
of ultracompact jets. We discuss possible applications of 
opacity effects to studying the conditions existing in cen- 
tral regions of active galactic nuclei. 

Key words: galaxies: jets - quasars: general - quasars: 
individual: Cygnus A, 3C 309.1, 3C345 



1. Introduction 

Extragalactic jets are believed to be formed in the vicinity 
of an accretion disk around a supermassive black hole in 
the center of active galactic nuclei (AGN). Pressure and 
density gradients, as well as a toroidally shaped magnetic 
field produced by the rotation of the disk, may coUimate 
the outflowing material into two bipolar streams, and ac- 
celerate it up to relativistic speeds (see reviews by Begel- 
man, Blandford, & Rees, 1984; Begelman 1995). Physical 
conditions in the immediate proximity of an accretion disk 
determine many aspects of jet evolution at larger spatial 
scales. The extent of this ultracompact fraction of the jet is 
estimated to be ~ 0.1 — 10 pc. Roland et al. (1994) model 
it as a turbulent region of < 1 pc in size. Marscher & Gear 
(1985) estimate the projected size of ultracompact jet to be 
~ 0.3/i~^ pc, based on VLBI observations of a radio flare in 
3C273. Konigl (1981) and Zensus, Cohen, & Unwin (1995) 
regard an unresolved core as an ultracompact jet. 

In images obtained with Very Long Baseline Interfer- 
ometry (VLBI) , the core is usually identified with the most 



Send offprint requests to: A. P. Lobanov 



compact (often unresolved) feature exhibiting a substantial 
flux and flat spectrum across the radio band. At any given 
frequency, the core is believed to be located in the region of 
the jet where the optical depth is r = 1. The core absolute 
position, rcorc, should therefore depend on the observing 
frequency, Uahs- Konigl (1981) gives rcore oc (i^obs)"^^''' ■ 
The power index kr depends on the shape of electron en- 
ergy spectrum and on the magnetic field and particle den- 
sity distributions in the ultracompact jet. Observed offsets 
of the core position at different frequencies have been re- 
ported for several sources including 1038-1-528 A (Marcaide 
& Shapiro 1984), 4C 39.25 (Guirado et al. 1995), 3C395 
(Lara et al. 1996), and 3C 309.1 (Aaron 1996). Several stud- 
ies of the core position offset have been undertaken for 
3C 345 (Biretta, Moore, & Cohen. 1986; Unwin et al. 1994; 
Zensus et al. 1995). 

In this paper, we discuss synchrotron self-absorption 
and free-free absorption in the nuclear regions of AGN. We 
use the frequency dependence of the VLBI core position as 
a tool for determiningthe physical conditions of ultracom- 
pact jets. In section Q, we describe a model adopted for 
ultracompact jets, and outline the relations between core 
shift and physical properties of the jets. Measurements of 
the shift of VLBI core in radio sources are discussed in sec- 
tion ^ In section ^, the measured core offsets are applied 
to deriving the magnetic field distribution and physical 
conditions in the central regions of Cygnus A, 3C 309.1, 
3C345, 3C395, 4C 39.25, and 1038-^528 A. 

Throughout the paper, we use a Hubble constant Ho = 
100 /ikms~^ Mpc~^ and deceleration parameter qo = 0.5. 
Unless defined otherwise, all quantities are in the cgs units, 
except for distances which are given in parsecs. 

2. Ultracompact jets 

A region of significant particle acceleration is assumed to 
exist at the very center of an AGN. The accelerated parti- 
cles are confined by an ambient plasma with steep pressure 
and density gradients along the rotational axis of the cen- 
tral engine. A highly coUimated relativistic outflow (jet) 
can be formed in result. The jet can contain both protons 
and electron-positron pairs (Sol et al. 1989). Interactions 
between the jet and the ambient medium are usually ig- 
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nored, and the jet itself is assumed to have no transverse 
velocity gradient. The jet hydrodynamics is then given 
by a stationary adiab atic flow satisfving Berno ulli's equa- 
tion (v • V)7jp^/* = |(Daly fc Marschor 1988)| . A nozzle is 



formed at a distance rc where the pressure p drops below 
4/9 of its initial value, and the flow becomes supersonic, 
with bulk Lorentz factor 7j(rc) = 1.22. Beyond the noz- 
zle, the flow is accelerated by the conversion of internal 
relativistic particle energy 7e to bulk kinetic energy 7kin, 
so that 7kin oc 7e oc r'^ (e =const). The resulting parti- 
cle energy distribution N{'yc) can be approximated by a 
power law: N{ye) = Noy'" for 7min(r) < 7e < 7max(r-), 
where s is the particle energy spectral index. The jet radi- 
ation is then described by the inhomogeneous synchrotron 
spectrum with spectral index a = (1 — s)/2. Unless other- 
wise specifled, the jet bulk Lorentz factor 7j, and viewing 
angle 6 are assumed constant. The jet geometry can be ap- 
proximated by a conical or paraboloidal expansion, with 
the transverse radius R — ar^ [e < 1), for rc < r < rmax, 
where rmax can be taken ~ 100-1000 rc, for most of parsec- 
scale jets. The jet opening angle = arctan(i?/r) is con- 
stant for e = 1. The magnetic field and particle density de- 
crease with r, and can be approximated as: B — B\{r\/r)"^ 
and A*' = Niiri/r)"^ , where B\, Ni are the magnetic field 
and the electron density at ri = 1 pc. Konigl (1981) shows 
that the combination m = 1, n = 2 can be used to de- 
scribe the observed X-ray and synchrotron emission from 
the most compact regions of VLBI jets. 

2.1. VLBI core and uHracompact jet 

In the scheme described above, the protons become subrel- 
ativistic in the rest frame of the flow outside a characteristic 
radius ro, and the bulk Lorentz factor freezes at the value 
~ 7j — 7(ro). This location in the jet is mos t likely to be 



From (|l|), fcr = ((3 - 2a)m + 2n - 2)/(5 - 2a) and 
fcb = (3 — 2a)/(5 — 2a). For the equipartition value of 
kr — 1, the choice oi m — 1, n — 2 appears to be the 
most reasonable. The corresponding fcr in this case does 
not depend on spectral index, and both B{r) and N{r) 
decline (taking m = 2, for instance, would result in n = 
0.5-|-a, which is unrealistic: N{r) would then remain nearly 
constant, or even grow, with increasing r). 

If the apparent core positions measured at two frequen- 
cies 1^1, i/2 (j^i < V2) differ by Armas, one can introduce a 
measure of core position offset 



D.r^ = 4.85 • 10" 



) Armas -Pl 
(1 + Z)2 



1/fcr l/fcr 



(4) 



where _Dl is the luminosity distance to the source, and 
frequencies are given in Hz. Qrv can be used for assess- 
ing the quality of core position measurements. For ideal 
measurements, firf = {B'^^ Fy^'^' sin 9 — const, for all fre- 
quency pairs. A dispersion in the values of fir^ derived from 
different frequency pairs reflects the inaccuracy of core po- 
sition data. If fcr has been determined inaccurately, or if 
it varies along the jet, then fin, will exhibit a systematic 
trend. The variations of fcr can occur when the jet crosses 
a region with rapidly changing absorption properties (e.g. 
the broad line region). In this case, fin, can be used to 
estimate the change of fcr. If two values Qrvi.2 (measured 
between frequencies vi and 1^2) and flrv2,3 (between U2 and 
1^3) are different, the relation between the corresponding 
fcr's is: 



kr2.3 ~ krl,2 [log (firi^l ,2 ) / log(r2r,yl ,3 )] 



(5) 



With a measured flrv, equation (|2|) can be applied to 
determine the magnetic field in the jet at r = 1 pc 



observed as the VLBI core (Marscher 1995) . It is also re- Bi = {fin / sin 



-jfcr/fcb^-l/fcb 



ferred to as the "injection point" (at which the relativistic 
plasma is assumed to be injected into the jet at 7j) in most 
of the models dealing with parsec-scale relativistic jets. 

At any given frequency, the VLBI core is observed at 
the location where the optical depth of synchrotron self- 
absorption is Ts = 1. For given A(r) and B{r), the corre- 
sponding Ta is (e.g. Rybicki & Lightman 1979): 



(6) 



and the absolute distance of the core from the central 
engine: 

-1 



rc. 



^{u) — fin, ^v^^'^' sin 6] 



(7) 



r{em + n — l) p^ + l 



(1) 



Here rric is the electron mass and e = 3/2 — a. The 
observed jet opening angle (j)o = (j)csc9. 6*2(0) is described 
in Blumenthal & Gould (1970), and C2(-0.5) = 8.4 ■ 10^°. 
Setting rs(r) to unity gives for the distance from the central 
engine to the core observed at frequency 1/ 



Formally, rcorc refers to the distance between the core 
and the sonic point, rc. However, since rcore ^ rc, in most 
cases, rcorc can be taken as a fair approximation for the dis- 
tance to central engine. It should also be stressed that rcorc 
in (^) depends only on measured quantities (Armas and fcr) 
and on the jet viewing angle 6 that can be determined from 
observations. Therefore rcore can be used as a reliable esti- 
mate of core position, or at least of the projected distance 
from the central engine (if 9 is poorly known). From (^ 
and (^), the magnetic field in the VLBI core observed at 
frequency v is given by 



r[pc] = {BI'^F/u 



with 



l/fcr 



F = (1 + z)-'[6.2 ■ 10'**C2(q)5;A-i,^o] 



(2) 



(3) 



BcoTc{v) ^ V 



l/fcr 



flriy 



(1 



"1/fcb 



(8) 



where C, = {k^/ki,) — m. For a typical a — —0.5, we have 
C = 0.5. 
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Table 1. Predicted core position shift between 5 
22 GHz 



Name 


Type 


Ar5,22 


L. 


yn 






[mas] 


[erg 




3C273 


LPQ 


0.78 


2.4 • 


10*^ 


3C216 


HPQ 


0.70 


3.1 ■ 


10« 


3C120 


LPQ 


0.61 


8.4- 


10« 


4C 39.25 


LPQ 


0.57 


2.2 ■ 


10** 


3C345 


HPQ 


0.37 


1.2 ■ 


10** 


1807+698 


BL 


0.26 


3.3- 


10*^ 


BL Lac 


BL 


0.11 


2.3- 


10*^ 




1 2 3 4 5 6 



2.2. Equipartition regime 

We now consider the equipartition between jet particle and 
magnetic field energy densities, with the magnetic field en- 
ergy density given 



Konigl 1979), where < 1, 



by kcAB'^/iT |(Blandford 

A — ln(7max/7min). In tMs case, kr — 1 (with m — 1, 
n = 2), and the core position offsets can also be used to 
determine the total (kinetic + magnetic field) power of the 
jet 



Lt 



2.03 • lO^'' [Qr 



,(1 + 2) 



3/2 e(3 + 2fc^A)2fA 



Sj sin 6 



(9) 



where 6 = ln(rmax/rc). The Lt 
derive the magnetic field at 1 pc 

1/4 



Bi « 2.92 • 10" 



can be further used to 



(10) 



The magnetic field in the core can be obtained, simi- 
larly to -Bcoro derived in (^). The equipartition can also be 
used for predicting the core offset in a source with known 
synchrotron luminosity Lsyn. In this case, the expected shift 
of the core position between two frequencies vi and 1^2 is 



Ci{l + z) U2 — 1^1 



f3{l ~ Pcos9)e 



2/3 



(11) 



with Cr = 4.56 • 10"^^, and assuming k^ = 1. Table |l| 
contains first order predictions of the core shift between 5 
and 22 GHz, for several prominent radio sources. 

The synchrotron luminosities in Table ^ are calculated 
from a database compiled by Ghisellini et al. (1992), using 
the model of Blandford & Konigl (1979) with 7max/7min = 
100, rmax/rc = 100. The jet opening angle (j) = 0.5° 
was used for all sources. The synchrotron self-Compton 
Doppler factors from the database are used for estimating 
the jet viewing angles and Lorentz factors, assuming the 
minimum jet kinetic power condition 7min = 1/sin^. This 
gives lower limits for synchrotron luminosities, and upper 
limits for the shifts. One can see, from Table 0, that the 
core shifts may be noticeably large, and can be measured 
from VLBI data. 



Fig. 1. Physical conditions in the jet in the presence of strong 
pressure gradients. Power index a describes the pressure gradient 
(p ex r"") along the jet axis. The corresponding magnetic field 
and particle density gradients have power indices mp (for m = 1) 
and n, respectively. The n and kr are given for different values 
of synchrotron spectral index a. 



2.3. External pressure gradients 

In the vicinity of an accretion disk, physical conditions in 
the jet become sensitive to the gradients of the pressure 
p oc r~" of the ambient medium. Since 7jP^^* — const, 
to satisfy Bernoulli's equation, the jet Lorentz factor and 
opening angle vary. Consequently, the corresponding gra- 
dients in B(r ) and N(r) should also change aloii g the jet 
in this region |(Georganopoulos fc Marscher 1996)| . The de- 
pendences of magnetic field and particle density distri- 
butions on pressure gradients are shown in Fig. |l[ If kr 
and spectral index a are measured, one can find, from 
Fig. hi the particle density and magnetic field distribu- 
tion satisfying the conditions of stationary adiabatic flow. 
The magnetic field at r ^ rc is oc r~'"p, nip = ma/4. 
The particle density distribution can be approximated by 
n = a(3- 2a)/4. 

If the jet is confined by gaseous clouds supported by 
thermal pressure from the central source and the stellar 
population with star density ps, the ambient medium den- 
sity can be modelled by an exponential decrease p(r) — 
po exp(— r^/d^), with the chara cteristic size of the cloud 
system d = {SkT/irGp^mpf -^ [(Blandford fc Rees 1974") 



Here mp is the protom mass, and T is the ambient plasma 
temperature. Such a distribution reflects the likely condi- 
tions in the broad line region (BLR) surrounding the nu- 
cleus of an AGN. The jet properties in this case are shown 
in Fig. |, for d = 400rc, a{rc) = 4, and m = 2. The pa- 
rameters are chosen so as to approach the equipartition 
at r ~ 3d. The jet opacity to synchrotron self-absorption 
increases significantly in the region of steepest pressure gra- 
dients. This implies that confinement effects may infiuence 
the observed frequency dependence of the VLBI core posi- 
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Fig. 2. Change of physical conditions along the jet axis, r. The 
pressure gradient a corresponds to gaseous clouds supported 
by thermal pressure and maintaining a mass distribution with 
spherically symmetrical gravitational potential. The cloud re- 
gion extends up to 400rc (rc refers to the distance at which the 
jet becomes supersonic). The equipartition regime is approached 
at the outer boundary of the cloud region, with mp = 1, n = 2, 
fcr = 1. Significant deviations from the equipartition are seen on 
smaller scales, resulting in stronger self-absorption in the inner 
parts of the jet (fcr < 2.5). 



tion in some radio sources. If ps and T are known, the core 
position measurements can be used for estimating the size 
of nuclear cloud system. 

2.4- Free-free absorption 

Foreground free-free absorption in a hydrogen plasma is 
expected to affect the VLBI-scale radio emission propa- 
gating through a dense nuclear environment. The existence 
of absorbing, circumnuclear plasma has been suggested in 
3C84 (Vermeulen et al. 1994; Walker et al. 1994), Cen A 
(Jones et al. 1996), and Cyg A (Krichbaum et al. 1997). In 
all three sources, only the counter-jet emission is believed 
to be absorbed. One can also expect to find absorption in 
the jet-side emission, if the jet viewing angle is relatively 
large, and the absorbing medium extends sufficiently high 
above the accretion disk plane. One possibility would be 
the broad line region f ormed by hydrogen plas ma clouds 
entrained from the disk 



(Cassidy 



Raine 1993)| . The opti- 
Levinson 



cal depth of free-free absorption is given by (e.g 
et al. 1995) 



rf(r) 



V gnc(r)ni{r)lpc , 



(12) 



where Ipc is the size of absorbing region. Assuming a 
pure hydrogen plasma with uniform density, we can take 
rie ~ Mi. The hydrogenic free-free Gaunt factor, g, can be 
evaluated numerically (Hummer 1988) . An analytical long- 
wave approximation of g given by Scheuer (1960) can also 
be used. 



The typical sizes of BLR inferred from reverberation 
mapping and modelling the optical broad line emission of 
AGN vary between 0.01 and 1 pc (e.g. Kaspi et al. 1996; 
Baldwin et al. 1996). The corresponding densities of hy- 
drogen plasma^_are_jif^~ 10* — lO^^cm"^. The disk-wind 
model 
up to ' 



(Cassidy & Raine 1993) predicts the BLR sizes of 
20 pc, with jiH between 10^ and 10^^ cm"^. There 



is evidence that weaker emission lines may originate from a 
very extended (r < 400 pc) thermal component with densi- 
ties riH ~ 1 0^ — 10* cm~^ and elec tron temperatures Tc ~ 



10* 



10" K 



(Ferguson et al. 1997). Most of these plasmas 



can alter significantly the GHz-range emission from VLBI 
jets (for instance, tj — 250 at 10 GHz, for = lO^cm"'^, 
T = 10" K, Ipc = 0.1). Udomprasert et al. (1997) report 
very high rotation measures, RM ~ 40 000radm~'^ in the 
VLBI core of the quasar OQ 172, which further supports 
the presence of high-density thermal medium around the 
ultracompact jets. 

For a spherical distribution of BLR clouds, we can take 
nc(r) cx enor"" (e is the volume filling factor of the cloud 
distribution). Then, from (p^, rt oc r"^""*"^. A crude es- 
timate n ~ 3 can be adopted (remembering that n refers 
essentially to plasma density variations along the jet axis). 
Then rcorc oc and the corresponding core shift is 

roughly 10 times smaller than that due to synchrotron self- 
absorption at typical VLBI observing frequencies. 



3. Measuring the core position offsets 

High-precision measurements of radio source absolute po- 
sitions which are required for alignment of VLBI im- 
ages at different frequencies are not readily available in 
most cases. Extensive absolute astrometry (see Fomalont 
1995) observations are necessary, in order to establish 
a reliable link between VLBI data at different frequen- 
cies. In quasi-simultaneous m ultifrequencv observations, 



the phase-referencing technique (Beasley fc Conway 1995) 



can be sufficient for the purpose of image alignment. If nei- 
ther of the abovementioned techniques is available, the fre- 
quency dependent shift of the core position can be deduced 
from comparison of observations made at close epochs, as- 
suming that the moving features observed in parsec-scale 
jets are optically thin and therefore should have their posi- 
tions unchanged. In this case, the offsets between the com- 
ponent locations measured at different frequencies will re- 
fiect the frequency dependent shift of the position of the 
source core. 

3.1. Core shift in 3C345 

An extensive long-term VLBI monitoring database is avail- 
able for 3C 345 (Zensus et al. 1995, 1997; Lobanov 1996). In 
the data for 3C 345, there are three close pairs of VLBI ob- 
servations: (1989.24 at 22.2 GHz)-(1989.26 at 10.6 GHz), 
(1992.44 at 22.2 GHz)-(1992.45 at 5 GHz), (1993.69 at 
5GHz)-(1993.72 at 22.2 GHz). In these pairs, the separa- 
tions between the observations in these pairs do not exceed 
10 days. We also use two other close pairs with separations 
of 51 days: (1992.71 at 8.4 GHz)-(1992.86 at 22.2 GHz), 



Table 2. Average offsets of the core position in 3C 345 



1 


2 


3 


4 


5 


Frequency pair 


Ar 


bAr 








[mas] 




[deg] 


[pc] 


5.0- 8.4 GHz 


0.05 ±0.03 


5% 


-70.8 ±4.1 


0.1 


10.6-22.2 GHz 


0.17 ±0.05 


20% 


-73.5 ±2.0 


0.6 


8.4-22.2 GHz 


0.21 ±0.06 


20% 


-71.4 ± 1.4 


0.8 


5.0-22.2 GHz 


0.33 ±0.10 


28% 


-77.1 ± 2.6 


1.2 



Notes: 2 - offset radius in mas; 3 - corresponding fraction 
of restoring beam; 4 - position angle of the offset direction; 
5 - projected linear distance corresponding to the offset. 



0.6 



a 0.4 



U 0.2 



0.0 



6 8 10 

Frequency [GHz] 
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and 70 days: (1993.69 at 5 GHz)-(1993.88 at 8.4 GHz). We 
measure the offsets in the closest and brightest components 
which have most reliably measured positions. For the pairs 
with 51 and 70 days separation, the positions of jet compo- 
nents have been corrected for the proper motions measured 
from the polynomial fits to the components trajectories 
(Lobanov 1996). 

The average angular offsets are compared in Table |^ 
with the restoring beam sizes of spectral index maps made 
from the corresponding frequency pairs. The direction of 
the core shift ( —75°) is similar to the position angle of the 
innermost part of the jet observed at 43 GHz (Krichbaum 
& Witzel 1992). Between 5 and 8.4 GHz, the position shift 
cannot be measured reliably because of the insufficient reso- 
lution and errors due to the proper motion corrections. The 
measured shift constitutes only about 5% of the restoring 
beam. With increasing frequency separation, the position 
shift becomes a prominent fraction of beamsize (up to 27% 
for spectral index maps between 5 and 22.2 GHz). Fig. |^ 
shows the measured offsets with respect to the reference fre- 
quency (22.2 GHz). The 5-8.4 GHz measurements are not 
included. The dashed line in Fig. H represents the best fit: 

^ rv ,/-1.04±0.16) 



3.2. Effect of the reference point offset on spectral imaging 



Whenever an offset constitutes a significant fraction of 
restoring beam, it can influence substantially the derived 
spectral properties. An example of this effect is shown in 
Fig. ^. The 5-22 GHz spectral index profiles shown in Fig. 
^ are measured along the jet ridge line in the nuclear region 
of 3C 345. The core extends from ±2 to ±5 mas. The profile 
obtained by aligning the images on their respective core po- 
sitions shows a strong gradient, with optically thin spectral 
index across a significant fraction of the core. Applying the 
measured core position offset between 5 and 22.2 GHz levels 
the spectral index across the entire core, with acoro ~ 0.2. 
Similar corrections applied to several other spectral index 
maps of 3C 345 have resulted in decreased peak values, and 



Fig. 3. Frequency dependent shift of the core position in 3C345. 
Dashed line corresponds to r oc i/"^""*. Reference frequency is 
22.2 GHz. 




Fig. 4. Spectral index profiles taken in the nuclear region of 
3C345 (5-22.2 GHz, September 93). Dashed line is the spectral 
index profile obtained by aligning both images by their respec- 
tive core positions. Solid line is the profile obtained after apply- 
ing the measured position offset of the core. 



smoother spectral index distributions across the nuclear re- 
gion. 

4. Properties of ultracompact jets 

To illustrate the implications of frequency dependence of 
VLBI core position, we calculate physical properties of ul- 
tracompact jets in 3C 345 and 5 other AGN with reported 
core position offsets. Table ^ summarizes the model param- 
eters (panel A), and gives the measured and derived quanti- 
ties (panel B). For all sources, we assume m = 1, n = 2, and 
use crude estimates 7max/7min = 100, and rmax/^c ~ 100 
(since only the logarithms of these ratios affect the cal- 
culated luminosities). The value of A'^i is determined in 
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Table 3. Model and derived parameters of ultracompact jets 



A. Model parameters 



Name 


z 


7j 


[deg] 


[deg] 


"j 


[cm^^] 


Refs. 


Cyg A 


0.0562 


5.0 


80.0 


7.0 


-0.5 


8000 


1,2 


3C309.1 


0.905 


5.0 


20.0 


2.0 


-0.6 


17000^ 


3 


3C345 


0.594 


5.6 


7.6 


2.4 


-0.6 


1300^ 


4,5 


3C395 


0.635 


7.1 


8.7 


0.5 


-0.6 


1800^ 


6 


4C39.25 


0.699 


11.0 


30.0 


4.0 


-0.5 


lOOOt 


7 


1038+528A 


0.678 


4.0 


20.0 


5.0 


-0.6 


1900^ 


8 



References: 1 - Krichbaum et al. 1997; 2 - Carilli & Barthel 1996; 3 - Kus 1993; 4 - Zensus et al. 1997; 5 - Lobanov 
1996; 6 - Lara et al. 1994; 7 - Alberdi et al. 1993; 8 - Marcaide et al. 1994. 



B. Measured and derived parameters^ 



Name 










' core 


itot 




Tb max 




[GHz] 


[pcGHz-i] 


[G] 


[G] 


[pc] 


[10^« ergs-i] 


[10*^ ergs-1] 


[10" K] 


Cyg A 


22/43 


2.1 ± 0.4 


0.07 ±0.01 


0.7 ±0.3 


0.10 ±0.02 


0.55 ±0.05 


0.07 ±0.01 


1.0 


3C309.1 


1.5/2.3 


16.3 ±7.2 


1.2 ±0.4 


0.6 ±0.4 


2.1 ± 0.9 


5.8 ± 1.7 


0.71 ±0.21 


7.2 




1.6/2.3 


15.4 ±8.8 


1.1 ±0.5 


0.5 ±0.4 


2.0 ± 1.2 


5.3 ±2.3 


0.66 ± 0.28 


7.2 




2.3/5.0 


7.2 ±3.6 


0.4 ± 0.1 


0.4 ±0.3 


0.9 ±0.5 


1.7±0.6 


0.21 ±0.07 


6.5 




5.0/8.4 


< 5.2 


< 0.2 


< 0.3 


< 0.7 


< 1.0 


< 0.1 


< 6.3 




8.4/15.1 


< 8.0 


< 0.4 


< 0.4 


< 1.0 


< 2.0 


< 0.2 


< 6.6 


3C345 


10.6/22.2 


13.1 ±3.8 


1.9 ±0.3 


0.4 ± 0.2 


4.4 ± 1.3 


19.8 ±3.2 


2.4 ±0.4 


25.8 




8.4/22.2 


10.8 ±3.1 


1.4 ± 0.6 


0.4 ± 0.2 


3.7± 1.0 


14.7 ± 2.2 


1.8 ±0.3 


25.1 




5.0/22.2 


8.1 ± 2.4 


0.9 ±0.2 


0.3 ±0.2 


2.8 ±0.8 


9.6 ± 1.6 


1.2 ±0.2 


24.2 


3C395 


2.3/8.4 


7.4 ± 2.4 


1.4 ± 0.6 


0.7 ±0.4 


2.2 ±0.7 


1.0 ±0.2 


0.12 ±0.02 


22.2 


4C39.25 


2.3/8.4 


6.3 ± 1.3 


0.6 ±0.2 


1.0 ±0.4 


0.6 ±0.1 


18.2 ± 1.6 


2.2 ±0.2 


2.4 


1038+528A 


2.3/8.4 


6.3 ± 1.2 


0.4 ± 0.1 


0.5 ±0.2 


0.8 ±0.2 


2.6 ±0.2 


0.32 ±0.03 


7.4 




8.4/15.1 


< 1.5 


< 0.13 


< 0.3 


< 0.2 


< 0.3 


< 0.04 


< 6.2 



Notes: ° - only the errors due to uncertainties of the core offset measurements are considered; f - derived by equating 
B\ to its equipartition value; % - Score and rcorc are given for v^hs ~ 22.2 GHz. 



most cases, by setting Bi to its equipartition value given by 
(px|). Synchrotron luminosity, Lsyn, and maximum bright- 
ness temperature, Tbmax, are calculated using the model of 
Blandford and Konigl (1979). For the purpose of compari- 
son between different sources, the offset measures, firi/ are 
given, rather than the angular offsets cited in the original 
publications. 



4-1- Cygnus A 

We adopt the component identification made by Krich- 
baum et al. (1997), and use their position measurements 
to derive the offset of the core in Cygnus A between 22.2 
and 43.2 GHz (observation epochs 1992.44 and 1990.40, re- 
spectively). Krichbaum et al. (1997) give Gaussian model 
fits describing the compact structure of the source. We use 
these models to calculate the shifts of the total of four 



features situated both in the jet (components J3, J4) and 
counter-jet (components CI, C2). The two closest jet-side 
components, Jl and J2, are blended together in the model 
at 22 GHz. The averaged shifts are: Arj = 74 /xas for J3 and 
J4, and Arcj = — 54^as for CI and C2. The negative sign 
for the shift on the counter-jet side indicates that the com- 
ponent separations from the core are shorter at 43.2 GHz. 
This fact further supports the identification, as well as fre- 
quency dependence, of the core location in Cygnus A. We 
then estimate that the shift of the core position amounts to 
64 ± 12 /xas, and use this value for calculating the physical 
properties of the jet. 

The estimated Bi — 0.07 G is consistent with the field 
stre ngth (-Bipc ~ 0.06 G, for (j>j = 7.0°) in a Poynting fiux 
jet (Lovelace & Romanova 1996) driven by a 3 • 10* M© 
black hole. Carilh & Barthel (1996, CR96 hereafter) use 
the unresolved core flux density of Cygnus A at 43 GHz, 
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and apply the minimum energy equation, to arrive at 
-Bm = 0.16 G in a nuclear region of 0.15 mas (0.11 pc) in 
size. We have r-43GHz ~ 0.05 pc, for the distance from the 
jet origin to the core observed at 43 GHz. With these val- 
ues, -Bi = 0.11-Bm tan(^j/r43GHz ~ 0.06 G. 

The total jet power, Ltot, can be related to the lumi- 
nosity of extended radio lobe, Lr: Ltot = Lr/*? (CR96), 
with ?7 =0.01-0.1 (Leahy 1991) describing the efficiency 
of converting bulk kinetic energy into radio luminosity. In 
Cygnus A, the lobe radio luminosity Lr = 4.4- 10*^ ergs~^ 
(CR96), which results in r; = 0.08, for our derived Ltot. 

The predicted location, rc.e., of the central engine in 
Cygnus A is offset by about 0.14 mas from the core. Com- 
parison between the 22 GHz flux densities of two bright fea- 
tures (J2 on the jet side, and CI on the counter-jet side) 
evenly separated from rc.c. gives 2.2 ± 0.3 for the jet-to- 
counter-jet ratio (this figure should be viewed as an upper 
limit, if a foregorund absorbing medium is present). The 
calculated ratio results in 6j — 80.6 ± 1.2 for the adopted 
7j = 5. Conversely, for 6 — 80.0 used in the model, the 
corresponding /3j = 0.92+°;?*. 
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Fig. 5. Opacity in the jet of 3C 309.1. Circles are the measured 
values of fcr in 3C 309.1 at different frequencies; solid line shows 
changes of fcr due to pressure gradients described in section 2.3. 
The model parameters are the same as in Fig. 2. 



4.2. 3C 309.1 

We use the position measurements of Aaron (1996 and priv. 
comm.), and adopt the jet geometry as modelled by Kus 
(1993). Only upper limits are available for the core offset at 
frequencies higher than 5 GHz, as the separation between 
the core and closest optically thin feature remains constant. 
The derived decrease substantially towards higher fre- 
quencies. We take it as evidence for larger fcr due to increas- 
ing opacity at shorter distances from the jet origin. Between 
1.5 and 2.3 GHz, kr is likely to be close to unity, so we take 
the corresponding derived quantities as best estimates of 
physical conditions in the jet. The estimated iVi = 17000 
is larger than in the other objects listed in Table pi possibly 
implying a stronger pressure confinement in 3C 309.1. The 
latter can be reconciled with the results from optical spec- 
troscopy (Forbes et al. 1990) suggesting that 3C 309.1 may 
have a massive [M > lOOOMoyr"^) cooling flow within a 
radius of ll.dh"^ kpc. 

The variations of kr necessary to account for the differ- 
ence in measured Qriy are compared in Fig. ^ with changes 
due to the pressure gradients described in section |2.3| . Here 
we assign rcoro(1.5 GHz) — 1200rc, treat upper limits of 
^Irv as measured points, and use the same model parame- 
ters as in Fig. ^ The behavior of jet opacity displayed in 
Fig. 1^ is in agreement with the measurements of Aaron et 
al. (1997), who report a moderately high rotation measure 
RM^oie ~ —1600 with strong gradients across a nuclear 
region of ~ 3h~^ pc in size. It should be stressed however 
that the model shown in Fig ^ is used for purely illustra- 
tive purposes, and we do not attempt to make quantitative 
conclusions or rule out other explanations for the observed 
core offsets in 3C 309.1. For instance, external free-free ab- 
sorption in the BLR clouds may result in similar changes 
of kr. 



4.3. 3C345 

Properties of the ultracompact jet in 3C 345 have been dis- 
cussed recently in several papers, based on X-ray (Unwin 
et al. 1994, 1997), VLBI (Zensus et al. 1995), and muhi- 
frequency observations of the source (Bregman et al. 1986; 
Webb et al. 1994; Stevens et al. 1996). We adopt a spectral 
index aj — —0.6 (Unwin et al. 1994), and derive the geo- 
metric parameters for our model from the statistical prop- 
erties of superluminal features embedded in the extended 
jet (Zensus et al. 1995; Lobanov 1996). 

To obtain a plausible estimate of Sj, one can assume 
that the jet carries the least kinetic power, so that 7j — 
(1 + Pa.pp)^'^ (the assumption is valid only if the observed 
speeds are those of bulk motions in the jet). Under this 
assumption, the variations of viewing angle are derived for 
the jet components in the immediate vicinity of the core 
(Lobanov 1996); and we estimate 7j « 5.6 and 9j ~ 7.6° for 
the Lorentz factor and viewing angle of the ultracompact 
jet in 3C 345. These values are close to the fitted values 
from the synchrotron self-Compton (SSC) models applied 
to 3C345 (Unwin et al. 1994; Webb et al. 1994). 

The jet opening angle, (jij « 2.4°, is calculated from 
measuring the jet size within 1 mas distance from the core. 
The particle density A'^i = 1300 cm~^ is almost 20 times 
smaller than the value given in Unwin et al. (1994) for an 
SSC model of the VLBI core. Zensus et al. (1995) have indi- 
cated that the SSC value of A'^i might be roughly 10 times 
smaller, to accommodate the magnetic field and particle 
density found in the jet superluminal features. The latter 
prediction is consistent with our estimate of A'^i within a 
factor of 2. 

The calculated offset measures in 3C 345 given in Ta- 
ble ^ decrease at lower frequencies. The corresponding kr 
must decrease by about 20%, between 5 and 11 GHz (which 
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Table 4. Derived parameters of the ultracompact jet 



V 


-^corc 


Arpc 




[GHz] 


[G] 






5.0 


0.09 ±0.05 


16.3 ±4.6 


0.57 ±0.16 


8.4 


0.15 ±0.08 


9.7 ±2.8 


0.34 ±0.10 


10.6 


0.18 ±0.10 


7.7 ±2.2 


0.27 ±0.08 


22.2 


0.38 ±0.21 


3.7 ± 1.0 


0.13 ±0.04 



Notes: the angular offsets are given for a jet with Q — 7.6 



is within the errors of the fit given in Fig. Assuming 
that A;r(5GHz) = 1, the decrease in Qrv corresponds to 
A:r(llGHz) = 0.81 ± 0.14. This implies that the optical 
depth in the jet also becomes smaller at higher frequen- 
cies — a situation that does not have a simple explanation 
in the frameworks of both synchrotron self-absorption and 
external free-free absorption. With this argument in mind, 
and also considering the magnitude of errors in the derived 
offset measures, it appears more likely that the decrease of 
Qrv results from blending effects due to the limited resolu- 
tion. With decreased angular resolution at lower frequen- 
cies, the blending between the core and a nearest jet feature 
(remaining unidentified in the VLBI maps) should become 
progressively stronger, resulting in a systematic trend in 
Qrv However, we cannot entirely exclude an explanation 
involving peculiar physical conditions in the source. We 
therefore take the average Qrv = 10.7±5.4 as the reference 
value for 3C 345. 

The corresponding average characteristic magnetic 
field, Bi — 1.4 ± 0.7 G, is lower than most of the values 
obtained by Webb et al. (1994), but can be reconciled with 
the results from Unwin et al. (1994), considering the dif- 
ference of the adopted particle densities. 

The average total power of the jet is Lt ~ (1.5 ± 0.4) • 
lO^'^ergs'S and we find Lsyn « (1.8 ± 0.5) ■ 10*^ergs-^ 

the maximum bright- 
10^^ K. The obtained 



With the derived magnetic field, 
ness temperature is Tbmax ~ 2.5 
value is consistent with the 10^^ K inverse Compton limit, 
similarly to the mean values of the maximum bright- 
ness temperature in a sam ple of superluminal sources 



[(Vermeulen fc Cohen 1994) . Arguments based on the hy- 
pothesis of equipartition between relativistic particle and 
magnetic energy density (ep/es ~ 1) predict a lower value 
of maximum brightness temperature Tbeq ~ 4.5 • lO'^^ K, 
Readhead 1994). The relativistic particle energy density, 
T^■^ and the energy density in the magnetic field, 
SB = By Sir oc T-\ This gives ep/es ~ (Tbmax/rboq)*''^ ~ 
2 ■ 10^ for the ultracompact jet, and therefore the jet 
must lose its energy through X-ray emission due to inverse 
Compton scattering. Unwin et al. (1994) argue that phys- 
ical conditions in the extended jet may as well deviate sig- 
nificantly from the equipartition. According to Unwin et 
al. (1994), Ep/es ~ 10^ - 10* at r ~ 10-100 pc, which can 
be a natural consequence of evolution and radiation losses 
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Fig. 6. Magnetic field distribution in 3C 345. Squares show the 
magnetic field in the compact jet derived from the frequency 
dependent shift of the core. Circles are the homogeneous syn- 
chrotron model estimates of magnetic field in the extended jet 
components (Zensus et al. 1995). Triangles show the character- 
istic magnetic field values from a model of magnetized accretion 
disk (Field and Rogers, 1993). Diamonds are the theoretical es- 
timates from Kardashev's (1995) model for the dipole magnetic 
field around a supermassive rotating black hole. 



in a plasma with original Sp/sB ~ 10^. 

4.4- Magnetic field m the jet of 3C 345 

Table ^ summarizes the derived positions and magnetic 
fields in the VLBI core of 3C 345 at different frequencies. 
According to the derived positions, the location of the core 
observed at 22 GHz should be about 4pc away from the jet 
origin. The angular offsets are all smaller than the VLBI 
beam size at corresponding frequency; therefore, the central 
engine is likely to be blended with the VLBI core even if the 
opacity is low enough for the emission from the innermost 
regions of the jet to be observed. Under such circumstances, 
evidence for emission from the inner jet can only be found 
in the spectral index or turnover frequency maps (Lobanov 
et al. 1997). 

In Fig. ^, the derived magnetic field is compared with 
the values obtained from synchrotron emission calculations 
for the extended jet (Zensus et al. 1995), as well as with 
model predictions for the central engine in 3C 345. The 
theoretical values of magnetic field in the close vicinity of 
central engine are calculated from a model of thin, magnet- 
ically driven accretion disk (Field and Rogers 1993), and 
from a model of super massive black hole surrounded by 

The accretion 



a strong magnetic field (Kardashev 1995) 
— ^3 



and B 



160 G 



disk model predicts B « 10 at 
at r « 24_Rg (_Rg — 2GMbh/c^ is the Schwarzschild ra- 
dius of the central black hole). The mass of central black 
hole can be related to Bi, so that Mbh ~ 2.7 • l(f BiMq\ 
and in 3C345: Mbh « 4 • IO'^Mq. Kardashev (1995) sug- 
gests that a strong dipole magnetic field may exist in the 
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vicinity of a supermassive black hole. The maximum field 
strength can be estimated from the equilibrium relations 
Mbh ~ Md ~ Mb, where Mbh, M^, Mb are masses of 
black hole, its disk, and its magnetic field respectively. 
Then B„ax = c^G-^'^M-j^i = 2.4 • lO"(Mbh/M0)-\ This 
gives for 3C 345: -Bmax ~ 6.3 • 10^ G. For a dipole mag- 
netic field, one obtains B oc r~^), so that B ~ 0.8 G at 
1 pc, slightly lower than our estimate in Table ^ To match 
the derived Bi, the mass of central black hole in Karda- 
shev's model must be Mbh = {cB^^^r^^'^)/{2V2G^^^) ^ 
7 ■ IO^B^/V^/^Mq « 8.3 • lO'' Mq. 

4.5. 3C395 and 4C39.25 

The core offsets are known from geodetic VLBI measure- 
ments between 2.3 and 8.4 GHz (Lara et al. 1996; Guirado 
et al. 1995). We use the models derived for the parsec- 
scale jets in these sources (Lara et al. 1994; Alberdi et al. 
1993), and adopt the equipartition values of A'^i. The re- 
sulting Qrv are similar to the 2.3-5 GHz offset measure in 
3C 309.1. Other derived parameters are also comparable 
to their counterparts in 3C345 and 3C 309.1. The derived 
synchrotron luminosity of 4C 39.25 is quite high, compared 
to Lsyn (10-90 GHz) = 5.5 • lO^^ergs"^ given in Bloom & 
Marscher (1991), even taking into account the spectral lim- 
its. 

Recent observations of 4C 39.25 at 22.2 and 43.2 GHz 
(Alberdi et al. 1997) indicate no detectable position shift 
in the source. From the measured fin, between 2.3 and 
8.4 GHz, the expected shift between 22.2 and 43.2 GHz is 
less than 0.05 mas. The latter value is below the image res- 
olution at both frequencies (« 0.3 and « 0.17mas). 

4.6. 1038+528 A 

Marcaide & Shapiro (1984) reported a 0.7 mas shift be- 
tween 2.3 and 8.4 GHz in 1038-1-528 A. Subsequently, Mar- 
caide et al. (1985) derived 0.5 < K < 1.4, based on VLBI 
data at 1.4, 2.3, 8.4, and 10.6 GHz. Marcaide, Elosegui & 
Shapiro (1994) applied a relativistic shock model by Gomez 
et al. (1993) to explain the observed shift. Rioja & Porcas 
(1997) have made position measurements at 2.3, 8.4, and 
15 GHz, and argued that at least a fraction of the observed 
offset can result from limited resolution and blending at 
lower frequencies. The 0^1,(2.3-8.4 GHz) calculated from 
the data of Rioja et al. (1997) is similar to the values ob- 
tained for 4C 39.25 and 3C395 (although the projection 
effects should be kept in mind, in view of the different view- 
ing angles in these sources). The upper limit given for the 
offset between 8.4 and 15 GHz corresponds to > 4.4, 
making it rather difficult to be explained by synchrotron 
self-absorption. Even for free-free absorption, it would re- 
quire very strong density gradients (n > 4) in the absorbing 
medium, in order to reproduce the given limit on kr. 

5. Opacity in the jets 

Assuming that the frequency term in (^) dominates core 
position offsets, the measured flrv can be used for compar- 
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Fig. 7. Relative changes of opacity in the ultracompact jets. 
Reference point is the 1.5—2.3 GHz offset in 3C 309.1, with fcr = 1 
and r = 1200rc assigned to it. For the rest of the data, the fcr 
and r are calculated from the corresponding Qrv measured in 
the respective source rest frames. Dashed hne is a model similar 
to that described in section 2.3, but with a{rc) = 6. 



ing the relative opacity in different sources. To make such a 
comparison, the offset measures and observing frequencies 
must be transformed to the respective source rest frames: 
Q'j.^ — f2rySinSj(l + 2); v' — 1/(1 + z). We then postulate 
that the largest Q.'^^ corresponds to fcr = 1, r = 1200rc, 
and use it as a reference point for calculating the fcr for 
the rest of the offset measurements. The results are plot- 
ted in Fig. ^. In 3G 395, the calculated fcr is very large 
(fcr > 17), suggesting that the jet viewing angle may be 
greater than the value cited by Lara et al. (1994). The ap- 
parent increase of fcr at shorter radial distances is consistent 
with the self-absorption scenario described in section 2.3 
(the measured 57^^ become smaller at higher frequencies 
at which the opacity is larger). The linear scale in Fig. ^ 
serves only as an illustration, since the intrinsic proper- 
ties of the sources (reflected by the B^^F term) may be 
different. Detailed studies of each individual object are re- 
quired for establishing a connection between the jet linear 
scale and opacity properties. A general prediction is that 
the position offset should be larger in flat -spectrum cores 
in which fcr ~ 1. When fcr increases, the offsets are ex- 
pected to become smaller, and the core spectrum should 
be inverted. 

A reasonably good agreement found in Fig. ^ between 
the derived and model fcr indicates that synchrotron self- 
absorption may be responsible for the observed proper- 
ties of optically thick emission from ultracompact jets. 
One can expect then to find similar physical conditions 
in the VLBI cores of different sources, if the measurements 
are done at the same rest frame frequency. To illustrate 
this, we calculate the magnetic fields and core distances at 
v' — 22 GHz, and plot them versus jet synchrotron lumi- 
nosity (Fig. g). The derived rcore increase, with increasing 
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Fig. 8. Magnetic field (filled circles) and distance of the core 
(squares) with respect to the jet synchrotron luminosities. The 
data are shown for v' = 22 GHz. Dot-dashed line 



Lsyn, and follow roughly the proportionality rcoro oc Lsyn 
resulting from (^) and The corresponding magnetic 
field however remains nearly constant, with an average of 
Bcorc = 0.30 ± 0.08 [G]. One may expect the values of 
magnetic field to vary stronger, and exhibit a rather large 
scatter, if jet emission is affected by external factors such 
as free-free absorption. The homogeneity of Bcoro seen in 
Fig. ^ suggests that the jet plasma is a primary factor de- 
termining the location and properties of VLBI cores, and 
that intrinsic physical conditions in the jets must be fairly 
similar. If this suggestion holds for other sources, the ex- 
pected magnetic field in VLBI core should be of the or- 
der of Bcoro(i^ob.) « 0.3 [0.045z/obB(l + [G], for r/obs 
measured in GHz. 



6. Conclusions 

The frequency dependent position shift of VLBI cores can 
be used for studying the most compact regions of extra- 
galactic jets, and obtaining quantitative estimates for phys- 
ical conditions in the immediate vicinity of central engine. 
In the absence of high-precision position measurements, a 
more simplistic approach can be employed, assuming that 
positions of optically thin features in the jet remain the 
same at all frequencies. Then, the position offsets of op- 
tically thin features can be interpreted as a frequency de- 
pendent shift of the self-absorbed core of the jet. In some 
sources (1038+528 A, 3C345, Cygnus A), insufficient res- 
olution and blending at lower frequencies may undermine 
the offset measurements. A way to detect and, to a cer- 
tain degree, to correct the errors due to the blending can 
be found in introduction of an offset measure, Q,rv, which 
is supposed to remain constant for the case of ideal mea- 
surements and unchanged opacity of the jet. Listed below 



are the main results from deriving the offset measures in 6 
sources with reported position offsets of VLBI core: 

1. The jet luminosity, the maximum brightness temper- 
ature, the core magnetic field and location with respect to 
the jet origin can be determined, from offset measures. For 
the well studied sources (Cygnus A, 3C345), the derived 
values are shown to agree well with values determined by 
other methods. In Cygnus A, the jet-to-counter-jet ratio 
determined from the estimated location of the central en- 
gine results in a self-consistent source geometry. In 3C 345, 
the derived maximum brightness temperature indicates 
that the ultracompact jet is strongly particle-dominated, 
and must release its energy through inverse Compton scat- 
tering. The obtained magnetic field distribution in the ul- 
tracompact jet of 3C 345 is in a good agreement with the 
values derived from inhomogeneous synchrotron models ap- 
plied to the extended jet components. Based on the de- 
rived characteristic magnetic field, it appears more likely 
that the jet is produced by a thin, magnetized accretion 
disk, rather than by a single supermassive black hole with 
a strong dipole magnetic field. 

2. External pressure and density gradients typical for 
the broad line region may change the optical depth along 
the jet via both synchrotron self-absorption and free-free 
absorption in an ambient medium. The offset measures de- 
rived for different frequency pairs and different objects can 
be used for studying the opacity effects in ultracompact 
jets. Changes of jet opacity in 3C 309.1 deduced from the 
variations of fcr appear to be consistent with a self-absorbed 
jet propagating through a region with strong density gradi- 
ents. Similar conclusion can be drawn from the variations 
of fcr determined for all studied objects, although the agree- 
ment with the model may in this case be coincidental (one 
argument for making such a reservation is that the initial 
pressure gradients required to fit the data for all sources 
simultaneously are very high: N{r) oc r~®, at r ~ re). 

3. The distance at which the observed VLBI core is lo- 
cated at a given frequency is scaled with jet luminosity, so 
that Tcorc oc Ls/n- At the same time the magnetic field in 
the core appears to be almost unchanged, from one source 
to another. This can be taken as another evidence for the 
synchrotron self-absorption to be a primary factor infiu- 
encing the observed properties of the jet core. 

Further, more detailed studies of absorption in ultra- 
compact jets are required for confirming the conclusions 
stated above, and for constructing better models of indi- 
vidual sources. Nearly simultaneous, multifrequency VLBI 
data obtained with phase-referencing should be most useful 
for a reliable determination of the core offsets, and subse- 
quent studies of the most compact regions in active galactic 
nuclei. 
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